Effects of extended impurity perturbation in d-wave superconductor 
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We describe the effects of electronic perturbation distributed on nearest neighbor sites to the 
impurity center in a planar d-wave superconductor, in approximation of circular Fermi surface. 
Alike the behavior previously reported for point-like perturbation and square Fermi surface, the 
quasiparticle density of states p(e) can display a resonance inside the gap (and very weak features 
from low symmetry representations of non-local perturbation) and asymptotically vanishes at e — > 
as p ~ e/ln 2 e. The local suppression of SC order parameter in this model is found to be somewhat 
weaker than for an equivalent point-like (non-magnetic) perturbation and much weaker than for a 
spin-dependent (extended) perturbation. 

PACS numbers: 74.25.Jb, 74.62.Dh, 74.72.-h 



INTRODUCTION 

The study of the density of states (DOS) in high-T c 
superconducting (SC) metal oxides has motivated many 
theorists and experimentalists through the last years, 
because it defines such fundamental physical parame- 
ters as the quasiparticle conductivity a, the penetration 
length A, the electronic specific heat, etc. This study is 
guided by the facts that i) the charge carriers are prac- 
tically confined to the Cu02 planes [l| and character- 
ized by 2D wavevectors k = (k x ,k y ) and ii) the SC or- 
der parameter has d-wave symmetry 2j with four nodal 
points ki = (±fcp/ v2, ±kp /v2) > where the Fermi surface 
crosses the nodal directions k x — ±k y , the SC gap func- 
tion Ak turns zero, and the quasiparticle dispersion law 
= v£w + Aw coincides with that of normal metal, £k- 

materials 



Another important factor is that the high-T c 
are the so-called doped metals, where the Fermi energy sp 
is defined by the density of charge carriers, introduced by 
the doping process Q . This very process creates the scat- 
tering centers for quasiparticles, due to random Coulomb 
fields from ionized dopants. Other scatterers, not related 
to the density of carriers, can be additionally introduced 
into the system, and all them can produce considerable 
effects on its physical properties |4j . 

In particular, some resonances can emerge in the quasi- 
particle spectrum with d-wave gap symmetry 0, @ , even 
at low concentrations of impurities. Such resonance man- 
ifests itself in a maximum of DOS at a certain energy 
e res < A = maxk |Ak|, as well as in logarithmic suppres- 
sion of DOS at £ -C £ T es (the energies being referred 
to e-p)- These conclusions can be directly compared to 
experimental results, as those obtained in the scanning 
tunnelling microscopy experiments 

The perturbation that impurities introduce into the 
electronic subsystem of crystal, depend either on their 
positions with respect to the lattice and on the potential 
they produce on nearest matrix sites. Within the sim- 
plest possible model, where an impu rity only disturbs a 
single site in the lattice ja, la, 111 [lCJ, llll LL2| , the poten- 



tial is characterized by a single perturbation parameter. 
This point-like perturbation model allows one to obtain 
a simple solution for the quasiparticle DOS in terms of 
their Green functions, leading to the above mentioned 
possibility of low energy resonances. However, in reality, 
the impurity perturbations in high-T c materials are not 
exactly point-like but rather extended to a finite number 
of nearest neighbor lattice sites to the impurity center. 
This raises an important question on how robust are the 
results of point-like approximation to the spatial extent 
and geometry of impurity perturbation. The opposite 
limit to the point-like perturbation, when the defect is 
much bigger of the Fermi wavelength and can be treated 
quasiclassically , hardly applies to real atomic substi- 
tutes in high-T c systems where the perturbation extends 
to few nearest neighbors of the impurity site. Recently, 
an example of such extended impurity center was consid- 
ered for a specific type of spin-dependent perturbation 
|l4j . Here we use a similar approach to get comparison 
of the point-like and extended perturbations within the 
same type of perturbation operator, impossible for the 
above mentioned case. 

Usually, the treatment of low energy excitations in d- 
wave systems involves a certain parametrization of the 
spectrum in the vicinity of nodal points. Thus, in the 
popular approach proposed by P. Lee , one expands the 
difference k — k; in the local axes, e^i = hi/kp and e^2 = 
(ei+i,i — e;_i : i)/2 (Fig. QJi), and then approximates the 
spectrum components as: 



£k = Hvpki, 



Hvoko 



with two characteristic velocities t>2 <C vp. However, this 
implies a "square" geometry of Fermi surface and, strictly 
speaking, it only holds close enough to the half-filling 
condition Ep/W sa 1/2 (where W is the bandwidth). In- 
stead, in many high-T c compounds we have ep/W -C 1 
and the geometry of Fermi surface is closer to circular 
(Fig. H 3 ). In this case, a more adequate parametriza- 
tion of spectrum is obtained with £k = fovp (k — fcp) and 
Ak = A cos 2(fk0(e'^ ) — where tpk — Bxctan(k y /k x ) and 
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FIG. 1: Fermi surfaces for a d-wave superconductor, a) 
" Square- like" geometry at closeness to half- filling (ef/W = 
0.475). The local axes e\ t \ and ei,2 are shown explicitly for 
the wave vector k near the nodal point ki , in the rest of quad- 
rants they are related to the respective nodal points in an 
analogous way. b) "Circular-like" geometry at lower doping 
(ef/W = 0.125), almost indistinguishable from exact circle 
(thin dashed line). 



the theta-function factor restricts SC pairing to the BCS 
shell of width £d (Debye energy) around ev- 

Most of the known treatments of impurity effects in 
doped and disordered d-wave SC systems, including a 
self-consistent T-matrix analysis of DOS, were developed 
within point- like perturbation model and "square" geom- 
etry of Fermi surface 0, Q • The present paper is aimed 
to extend those studies, in order to include a more real- 
istic features, either of the impurity perturbation (which 
can affect several equivalent neighbor sites to the impu- 
rity ion) and of the Fermi surface geometry (in a more 
adequate circular approximation). 



PHYSICAL DESCRIPTION OF THE SYSTEM 



We use the Nambu spinors \l/ k = (et k a_ k x) where 



and ak,(T are the Fermi operators for quasiparticles 
with wave vector k and spin a and the model Hamilto- 
nian for disordered d-wave superconductor: 



(1) 



H — Hq + Hi mp , 

ff o = E*k(ekr 3 -A k n)* k , 

k 

H ™p = ~Jf £ e^'-^P^e^'-^^yv^. 

k,k',p S 

In what follows we denote by hats the matrices in Nambu 
indices, e.g., the Pauli matrices % and the matrix V — 
Vi mp T3 which describes the quasiparticle scattering by ex- 
tended (attractive) perturbation Vi„ lp around an impu- 
rity center p, over its near neighbors 5 (Fig. Formally, 
this perturbation only differs by the presence of T3 factor 




FIG. 2: Extended perturbation over four nearest neighbor 
sites to the impurity ion (its projection onto the CUO2 plane 
is shown by the dashed circle at the origin). 



from that considered in Ref. ^4|. The concentration of 
randomly distributed centers c = iV _1 ^ p 1 (where N is 
the number of cells) is supposed small, c -C 1. 
We define the Green function (GF) matrices as: 



G k ,k'(e) =«*k|4<» 



5 i(-<o)*(|$ k ( t ),$t #(0) j )dt) (2) 



where (. . .) is the quantum-statistical average with 
Hamiltonian, Eq. ^ and {a(t),b(0)} the anticommuta- 
tor of Heisenberg operators. The energy e is referred to 
the chemical potenial /i (identified in what follows with 
the Fermi energy ep). The observable characteristics fol- 
low from the averages of products of these operators (at 
given inverse temperature 0) , expressed through the cor- 
responding GF's by the spectral theorem: 

^ = r i8FZTMw a >>- (3) 

J — 00 ^ ' 

For the disordered system with the Hamiltonian, Eq. ^ 
we calculate GF's from the basic equation of motion: 

Gk,k' = G k <5k,k' — 

1 J- S^'^a^a^GiVG^y. (4) 



k",p,j 



Here G° 



(e — ^k^3 — AkTi) , and we expanded the 
structural function for impurity scattering in Eq. ^ as: 
e i(k'-k)S = y^ =1 a ik a ik ,. The functions 



cik x aky 
ai,k = 2 cos — ^— cos — 



a 3:k = 2 sin — cos 



a 2 ,k = 2 cos — sin 
a 4 ,k = 2 sin — suiy, 



realize irreducible representations of the C4 point group 
(j = 1 being related to A-, j = 2, 3 to E-, and j = 4 to 



3 



-B-representations [T^) and thus satisfy the orthogonality 
condition 



k - Sjji . 



(•5) 



The impurity effects on quasiparticlc spectrum are then 
naturally classified along these representations, alike the 
known results for magnetic impurities in ferro- and anti- 
ferromagnetic crystals 0, . 



T-MATRIX SOLUTIONS FOR THE GREEN 
FUNCTIONS 

The orthogonality of the a^k functions implies that in 
the general solution to Eq. 0] 



G k = 



-Si 



(6) 



the self-energy matrix turns additive in these represen- 
tations: Ek = XI j ^jk- Each partial term in the latter 
sum can be given by a specific group expansion (GE), 
like those known for point-like impurity perturbations in 
normal or superconducting 0, 0] systems and also 
for extended perturbations in magnetic systems |l7|. 



Sj k = —dTj 1 1 — cA 



cAj 



n#0 



+c > ( (n) e- lkn + A 4 ? (n) ) x 



x l-l 2 (n) 



(7) 



Here Tj — V(l + VGj) 1 is the (renormalized) partial 
T-matrix, and the matrices Aj (n) represent indirect in- 
teractions (in jth symmetry channel) between scattcrcrs 
at sites and n: 



Aj(n) = -Gj(n)Tj, G 3 {n) 



Y^„ikn 2 P< 



Gj = Gj (0) 



The sum ^ n ^o m El describes all the processes in- 
volving pairs of impurities, it implies averaging in random 
impurity configurations and hence runs over all the lat- 
tice sites n, the omitted terms are for triples and more of 
impurities. This defines a generalization of the GE ap- 
proach for extended impurity centers in superconductors. 
If the series in the brackets is restricted to its first term, 
the self-energy matrix £k becomes independent of k: 



>£5 



(8) 



Then, for small enough concentration of impurities, the 
renormalization of T-matrices can be neglected and we 



arrive at 



Tj if = vii + cyy 1 , 



3 TV ^ 
k 



r<0 
j,k b k- 



(9) 



The matrix functions G° can be expanded in the basis of 
Pauli matrices 



G" = po (gjo + 9jin - gj 3 r 3 ) 



(10) 



where po = 4/(?rW) is the constant DOS in 2D nor- 
mal system (the absence of T2 component in Eq. ^| is 
related to the fact that the gap function Ak is chosen 
real). Passing in Eq. Elfrom summation in k to integra- 
tion in "polar" coordinates £k = £ and v?k = <p (Fig. lb) 
accordingly to the rule: 



_ V f ~ — 



-in 



2/po—H r2-n 
-u JO 



we calculate the dimensionless coefficient functions gji 
Some of them are zero by the symmetry reasons: gn = 
<74i = 0, the rest can be approximated as: 



g 3 o - ^ 

321 



a] go, 



-931 



9j3 ~ 



(11) 



Here a 2 are the average values of a 2 k over the Fermi 



surface: 



4(l-w), 



*2,3 



4u, 



2cj 2 , where 



the band occupation parameter uj = fi/W is supposed 
small, in concordance with the chosen circular geometry. 
The functions go and gi are known from the studies of 
point-like perturbations |20l l2lj| : 



i-2/po-fM i-2-k 



dtp 



£ 2 - A 2 cos 2 2<p 
e [l/p, - Fi (1 - £ 2 /A 2 ) /A] , (12) 



gi (e) 



A 
27 



cos 2 2<p dip 



o ./o e 2 - £ 2 - A 2 cos 2 2</j 
« A/e D + 2 [F 2 (1 - e 2 /A 2 ) 

+ (e 2 /A 2 )F!(l-e 2 /A 2 )] , (13) 

they include /i = /i(l — 2lu/tt) » /i and the functions 
Fi(z) = K(l/z)/^andF 2 (z) = Vz E (l/ z ) with full ellip- 
tic integrals of 1st and 2nd kind K and E, having similar 
behavior with the elementary functions obtained within 
square geometry |]j . In the same similarity, the function 
.93 = (47T)- 1 f_>^ tdt ft* dp/(e* - e - A 2 cos 2 2p) is 
practically constant: 173 w In ^/tt/(2lu) — 1, within the 
relevant energy range |e| <^ W, p. 

Using these results, we readily calculate the partial T- 
matrices, Eq. The most important contribution to S 
comes from the j — 1 term (^4-representation): 



1 i - 



VA VAgo - T 3 



a\po 



D 



(14) 



4 



where v A = a\V lmp po/(l - alV lmp p g 3 ) is the dimen- 
sionless perturbation parameter in the yl-channel, and 
Da(e) = 1 — v\gQ(e) is the energy dependent denomina- 
tor. In particular, it can produce a low energy resonance 
at e = e res such that ReDA{e r es) — 0, analogous to 
the above mentioned resonance from point-like impurity 
center. This requires that va exceeds some critical value 
« 2/tt. 

The contributions from j = 2, 3 (^-representation) are: 



'fO _ "I 
1 2,3 — — 



ve v E (go T 9in) - % 



D e 



(15) 



with the respective perturbation parameter ve = 
a\Vim P po/{l - atVimppags) and denominator De = 1 — 
v e \9o ~ Si) ■ ^ is l ess probable to have a resonance ef- 
fect in this channel at low occupation w < 1, since i) 
the parameter ve is reduced vs the j4 -channel value, and 
ii) there is a competition between Rc^q and Reg^ in the 
denominator De- 

The 5-channel contribution (j = 4) has the same 
structure as the ^4-channel term, Eq. ^] but with va re- 
placed by a strongly reduced value vb — a|^ mp po/(l — 



am 



impP0d3 ) 



4 v imph 

hence it turns even less important than the 



i?-channel terms. 
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FIG. 3: Density of states in the <i-wave superconductor with 
extended impurity centers (the solid line), for the choice of 
parameters W = 2 eV, p = 0.3 eV, e D = 0.15 eV, V imp = 0.2 
eV, c = 0.15. The arrow indicates the low-energy resonance 
by the yl-channel impurity effect and the dashed line repre- 
sents the pure d-wave DOS. Insets: weak i?-channel "antires- 
onances" at high energies (upper panel) and near the gap edge 
(lower panel). 



PERTURBATION OF OBSERVABLE VALUES 



Now we are in a position to describe the perturba- 
tion of basic observable characteristics of SC system 
by extended impurity centers. Thus, the global DOS 
is defined by the momentum diagonal GF's, p(e) — 
(TrAO^EklmTrGk, and, using Eqs. IgllHlHSlin Eq. 
H3 it is obtained as 

p(e) = ^-Img (e-Eo). (16) 

Here the scalar self-energy 

Pa \a\DA a\D E a\D B ) 

includes the effects of extended impurity centers in all 
three channels. The result of direct calculation from Eq. 
Hoi with use of Eq. El f° r the characteristic choice of 
parameters, W — 2 eV, p — 0.3 eV, £d = 0.15 eV, 
Vimp — 0.2 eV (giving for particular channels: va ~ 
0.934, v E ~ 0.088, and v B « 0.006), and c = 0.15, is 
shown in Fig. [3] It is quite similar to the known results 
for point-like impurities 0, 0, showing a reduction of 
the sharp coherence peak at e — A and emergence of 
a relatively broad low-energy resonance at e res (shown 
by the arrow), mainly due to the ^4 -channel effect. But, 
additionally, there are small "antiresonance" effects from 



the i?-channel (insets to Fig. at e ph A and at some 
high enough energy (~ 70 A in this case). Clearly, these 
i?-channel features shouldn't have any practical effect on 
the system thermodynamics. 

The local density of states (LDOS) on nth 
site is expressed in terms of GF's as p n (s) — 

(7rA0 _1 Ek.k' Im Tr e i(k ~ k,) ' n G k , k < and its variation 
^Pn(e) = Pn(s) — p(z)> compared to the mean value 
p(e) = N~ 1 ^2 n p Ii (e) (identical to the global DOS), is 
only given by the momentum- nondiagonal GF's: 

«p»(e) = 4? V ImTre*( k - k '>- n G k , k ,. (18) 

These functions are easily calculated for the simplest 
case of a single impurity center at p = 0: 

= ^J2 ",k^ k /:'^k^,k- , (19) 

3 

describing a finite effect on the local characteristics near 
the impurity. Thus, the quantity Sp n attains its maxi- 
mum value at n = 5, the nearest neighbor sites to the 
impurity. Using Eq. 1191 and the orthogonality relations, 
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FIG. 4: Local density of states on the nearest neighbor site to 
an extended impurity center, for the same choice of param- 
eters as in Fig. [3] (but supposing c — > 0). Note an overall 
enhancement of electronic density compared to remote sites 
from impurity (dashed line) and a much stronger effect of the 
low-energy resonance (the arrow). 



we expand this value in a sum: 
1 



-ik'-d 



k.k'j 



and present the overall maximum LDOS as: 
2po 



Pn=5(e) 



7T 

v A N A 



ol\D a 



Im \g (e)(l+ 

^ veN e vbN b \ 



a\D E 



\D B ) 



(20) 



Alike Eq. [T7|for global DOS and the case of Ref. pj for 
LDOS, the resonance contribution to Eq. 1201 comes from 
the j4-channel with the numerator N A = 2<?3 + v A {g$ + 
g 2 ), while other channels with Ne = 2g3+VE(go — <?? — <?f) 
and Nb = 2g 3 + VB{g$ + <jf) mainly contribute to renor- 
malization of the pure rf-wave DOS Pd{s) = 2/7rImgo(£)- 
The calculated from Eq. l2Ul behavior of LDOS on nearest 
neighbor sites to the impurity is shown by solid line in 
Fig. 2| It displays a low energy resonance (the arrow), 
much more pronounced than that in the global DOS, Fig. 
03 and an overall enhancement compared to the LDOS 
curve for remote sites from impurity p n ^oc — Pd (the 
dashed line). This picture can be compared with the 



direct experimental measurements of differential conduc- 
tance through the STM tip positioned close to and far 
from an impurity center 

In a quite similar manner, the local perturbation of 
SC order parameter can be considered. The local d- 
wave SC order in the unit cell containing the impurity 
(Fig . |3) is given by the average A 32 = 2V{as 3 ^as 2 ^) 
|14j. where V is the SC coupling constant and site oper- 
ators a njCr are expressed through band operators: a n , CT — 
^y-r/2 e lk n ak, CT - For Vi mp = 0, this average coincides 
with the uniform gap parameter: 



. 2V v-^ 
A = — > i 



>4 k-(<5 2 -<5 3 ) 



4AA r e v/^ 
— / [Fafl+xVFrtl + x 2 )]^, (21) 
71 Jo 

the latter expression (with the dimensionless d-w&ve cou- 
pling constant A = VpoLu) being obtained from Eqs. |3| 
and 1131 The integral in Eq. [2] behaves as a logarithm: 
f£[F 2 (l+x 2 )-F 1 (l+x 2 )}dx « (tt/4) ln(2.428a) at a > 1, 
thus providing A ss 2.428 £d exp(— 1/A). 

But for Vj mp 7^ this value is locally suppressed. The 
suppression is characterized by the dimensionless param- 
eter r] SU p = 1 — A32/A, confined between (pure SC) 
and 1 (complete local suppression of SC order) 0] , and 
it only results from non-diagonal GF's: 



N 



0-k,| a k,T/ 



sup 



2V 
AA 



E («-k,i«k', T )e l(k - fe - k '- 53) 



k.k'^k 



V 
4ttA 



,0 

J —' 



delm TrG"T°G"fi. (22) 



Using here Eqs. 1191 and ITU1 leads to the expression: 



4A r D/A 



1]sup 



71" JO 



F(x)dx, 



(23) 



where only the ^-channel terms contribute to the dimen- 
sionless function 

F(x) = --lm De 

(compared in Fig. [SJwith the integrand for the uniform 
gap equation, Eq. I21|) . Numeric analysis of this expres- 
sion for the above chosen perturbation parameters results 
in Vsup ~ 0.47. This is somewhat smaller than the re- 
spective value for point-like impurity in d-wave system 
6]: rj sup = l/(l + v 2 ) (assuming v equal to v A ), and more 
than twice weaker than almost complete suppression in 
the case of spin- dependent perturbation of the same di- 
mensionless magnitude jl4j . Those relations confirm the 
general Abrikosov-Gor'kov conclusion on the pair- 
breaking effects by impurities in superconductors, irre- 
spectively of their spatial extension. 



6 




e/A 



including the self-consistent ^4-channel T-matrix: 

f x (sc) = v[i + d[ sc) v]- 1 . 

Alike Eq. ^3 the self-consistent GF matrix = 
N^ 1 ^ k c^l^G^^ can be parametrized in Pauli matrices: 

G\ = PoQ-l (g-g 3 T 3 ). 

Then, using Eq. El hi Eq. |21 we readily conclude that 
j]( sc ) is diagonal in Nambu indices that is, within the 
considered self-consistent T-matrix approximation for ex- 
tended impurity centers, the scattering by dopants does 
not influence the <i-wave order parameter, the same as 
for point-like centers Q|. Hence, the self-consistency 
should be achieved only for the scalar function g — 

y(sc) 



(2N)- 1 Ek TrG k through the equation: 



7(e) = g {e 



00), 



(25) 



where the scalar self-consistent self-energy: 



FIG. 5: The dimensionless function F(e/A) (solid line) used 



in Eq. 



| to calculate the suppression parameter g S u P , at 



the same choice of parameters as in Fig. |1J compared to the 
integrand in the uniform gap equation, Eq. 1211 (dashed line) 
and its asymptotics 7rA/(4e) (dash-dotted line). 



SELF-CONSISTENT GENERALIZATION 

The above given analysis corresponds to the simplest 
restriction of the group series, Eq. [7| for self-energy to 
its first, single- impurity term with use of unperturbed 
GF's. The resulting linear approximation in impurity 
concentration: £ = — cY]j T®, is only justified when this 
concentration is low enough c < Co, where cq ~ pqSq is 
related to the characteristic energy scale £o for impurity 
perturbation (in this case £o ~ £res)- At higher con- 
centrations, c > Co, when perturbations from different 
impurity centers effectively overlap, a simplest way to 
take account of these collective impurity effects is pro- 
vided by replacement of £ by its self-consistent ana- 
logue £( sc ) = — cJ2jTj SC \ in the spirit of well-known 
self-consistent T-matrix approximation EJ, H3 • It was 
shown for the case of point-like impurity perturbation 
0, that effects of such self-consistency are most essen- 
tial at the lowest excitation energies, e <C A. In view 
of the similarity in the system response to point-like and 
extended perturbations and of the predominant role of 
the ^-channel at low energies, we can restrict the self- 
consistency of T-matrix only to its ^4-channel term. Then 
Eqs. I© and © are modified to: 



G 



(sc) 



(sc) 



(24) 



Po 



(sc) 



a\D E 



a\D B 



(26) 



includes the self-consistent denominator: D^ c ' (e) = 1 — 
v\g 2 (e). Then, passing to dimensionless energy x = e/A 
and denoting 



,(sc] 



°{g) = ag 



1 



2 2 

a g 



■0 



(27) 



with_a = v A , = cv A /{a 2 lPo A), (3' = c[2v 2 E /(a 2 D E ) + 
v B /(alD B )]/(v A p A) and 7 = A/[p(l - 2w/tt)], we ar- 
rive at the self-consistency equation for g — g(x) as: 



(x - a) {7- 



1 



1 - (at - try 



1 - {x - ay 



, (28) 



with a = a(g) defined by Eq. 1271 This equation is quite 
similar to that reported for point-like impurity centers 
and square- like geometry 0, , differing only by the ap- 
pearance of additional term /?' in Eq. |23 In the same 
way, Eq. 1281 admits two types of solutions in the en- 
ergy range of principal interest x — > 0. One of them, 
g = g^ ( x )j tends in this limit to a finite imaginary value: 
g^ x \x — > 0) — ► 170, defined by the equation 



-&f (7o) { 7 



xK 



y/1 + « 2 7o 7 2 (7o) 



1 + a^lf 2 (70) 



(29) 



with / (70) = 0' + 0/(1 + a 2 7o)- Another solution, g = 
g^ 2 \x), is vanishing in this limit: g^ 2 \x — > 0) — > 0, so 



7 



that all the denominators D in Eq. can be safely put 
equal unity, simplifying Eq. 1281 to: 



g = (x - a'g) I 7 - K 



1 — (x — a'g) 2 



(30) 



with a' — a(P+P'). Its solution has the same logarithmic 
asymptotics at x — > 0: 



g(x) w — 
a' 



1 



1 



7ra' ln(2i7ra'/ir) 



(31) 



as found in Refs. (J,|7j- Also the conclusion on the va- 
lidity range for each of the two solutio ns: far enough 
from and g^ close to the Fermi level |20j, remains true 
in the present situation. Thus the DOS p(e) at e <C e res 
should be suppressed as: 



2ttA [a' ln(27ra'A/e) 



(32) 



compared to its linear asymptotics p ~ p^ej A resulting 
from the linear in c approximation, seen in Fig. [3] 



CONCLUSION 

The Green function analysis is developed for the quasi- 
particle spectrum in a planar rf-wave superconductor 
with finite concentration of impurity centers which per- 
turb atomic energy levels on nearest neighbor lattice 
sites. A generalization of the method of group expan- 
sions for quasiparticle self-energy is obtained for such ex- 
tended impurity centers. The general picture of spec- 
trum restructuring is found quite similar to that pre- 
viously established for point-like impurity perturbation, 
though some new specific features due to extended nature 
of the perturbation are also indicated. In particular, it is 
found that the effects on the quasiparticle DOS and on 
the SC order parameter result from different irreducible 
representations of the point symmetry group of the im- 
purity center. The self-consistent procedure is proposed 
for higher concentration of extended impurities, general- 
izing the known formulations for point-like centers, and 
a qualitative similarity with that case is demonstrated. 
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